Application of the χ principle and unbiased predictive risk estimator for determining the regularization parameter in 3D focusing gravity inversion
نویسندگان
چکیده
The χ principle and the unbiased predictive risk estimator are used to determine optimal regularization parameters in the context of 3D focusing gravity inversion with the minimum support stabilizer. At each iteration of the focusing inversion the minimum support stabilizer is determined and then the fidelity term is updated using the standard form transformation. Solution of the resulting Tikhonov functional is found efficiently using the singular value decomposition of the transformed model matrix, which also provides for efficient determination of the updated regularization parameter each step. Experimental 3D simulations using synthetic data of a dipping dike and a cube anomaly demonstrate that both parameter estimation techniques outperform the Morozov discrepancy principle for determining the regularization parameter. Smaller relative errors of the reconstructed models are obtained with fewer iterations. Data acquired over the Gotvand dam site in the south-west of Iran are used to validate use of the methods for inversion of practical data and provide good estimates of anomalous structures within the subsurface.
منابع مشابه
Regularization parameter estimation for underdetermined problems by the χ principle with application to 2D focusing gravity inversion
Abstract. The χ-principle generalizes the Morozov discrepancy principle to the augmented residual of the Tikhonov regularized least squares problem. For weighting of the data fidelity by a known Gaussian noise distribution on the measured data and, when the stabilizing, or regularization, term is considered to be weighted by unknown inverse covariance information on the model parameters, the mi...
متن کاملAutomatic estimation of regularization parameter by active constraint balancing method for 3D inversion of gravity data
Gravity data inversion is one of the important steps in the interpretation of practical gravity data. The inversion result can be obtained by minimization of the Tikhonov objective function. The determination of an optimal regularization parameter is highly important in the gravity data inversion. In this work, an attempt was made to use the active constrain balancing (ACB) method to select the...
متن کامل3-D Projected L1 inversion of gravity data using truncated unbiased predictive risk estimator for regularization parameter estimation
متن کامل
3D Inversion of Magnetic Data through Wavelet based Regularization Method
This study deals with the 3D recovering of magnetic susceptibility model by incorporating the sparsity-based constraints in the inversion algorithm. For this purpose, the area under prospect was divided into a large number of rectangular prisms in a mesh with unknown susceptibilities. Tikhonov cost functions with two sparsity functions were used to recover the smooth parts as well as the sharp ...
متن کاملOn Concomitants of Order Statistics from Farlie-Gumbel-Morgenstern Bivariate Lomax Distribution and its Application in Estimation
‎In this paper‎, ‎we have dealt with the distribution theory of concomitants of order statistics arising from Farlie-Gumbel-Morgenstern bivariate Lomax distribution‎. ‎We have discussed the estimation of the parameters associated with the distribution of the variable Y of primary interest‎, ‎based on the ranked set sample defined by ordering the marginal observations...
متن کامل